Preparation of entangled states of two photons in several spatial modes 
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We describe a protocol capable of preparing an arbitrary state of two photons in several spatial 
modes using pairs of photons generated by spontaneous parametric down-conversion, linear optical 
elements and single-photon detectors or post-selection. The protocol involves unitary and non- 
unitary transformations realizable by beam splitters and phase shifters. Non-unitary transformations 
are implemented by attenuation filters. The protocol contains several optimization capabilities with 
the goal of improving overall probability of its success. We also show how entangled two-photon 
states required for quantum computing with linear optics can be prepared using a very simple and 
feasible scheme. 

PACS numbers: 03.67.Bg, 42.50.Dv 



I. INTRODUCTION 

Quantum information processing (QIP) has demon- 
strated an important development in recent years [ij. 
The significant advances in implementing various pro- 
tocols for QIP can, in the future, lead to long-distance 
quantum communication [2] or creation of a quantum 
computer Q. Among many experimentally tested phys- 
ical platforms, optical implementations appear particu- 
larly promising. Light is an ideal carrier of quantum 
information, as it can be transmitted over long distances 
through optical fibers or even in a free space. A crucial 
prerequisite for optical quantum information processing 
is the ability to generate and manipulate various highly 
non-classical and entangled states of light beams. Entan- 
gled states of light are necessary for instance for quantum 
teleportation [4, 5, 6, 7], quantu m g ates in quantum com- 
puting with linear optics @, 0, [l0|, [HI, , or quantum 
repeaters 

These wide potential applications motivate the ef- 
fort aimed at developing and demonstrating prepara- 
tion protocols for various interesting quantum states of 
light. Single and two-photon Fock states were condi- 
tionally prepared from two-mode squeezed vacuum by 
measuriiig the number of photons in the idler mode 
[13, m, E, [I3|- Schrodinger cat-like states were gen- 
erated by subtracting a single photon from pulsed or 
continuous squeezed beams [H, [l^ [13, [U, . Photon- 
added coherent states were produced experimentally us- 
ing down-conversion seeded with a weak coherent signal 
and conditioning on detection of a photon in the idler 
mode (23} . Schemes for preparation of arbitrary single 
mode states via repeated addition or subtraction of sin- 
gle photons have been suggested [H, [2^. A universal 
procedure for conditional preparation of arbitrary mul- 
timode states of light with linear optics, single photons, 
coherent states and single-photon detectors has been de- 
veloped [26j . 

A lot of attention has been paid to creation of two- 
mode N-photon entangled states, the so-called NOON 
states, [23, [M [H, [M [m, [13, m, m. it was shown 
theoretically and demonstrated experimentally that such 



states can be conditionally generated from input single- 
photon states by multiphoton interference in a properly 
tailored interferometer. In a similar way, projection onto 
arbitrary NOON- type state can be accomplished [ssj . 
The N-photon entangled states can find applications not 
only in quantum information processing, but can be also 
used for ultra-precise measurements [H, [13, [H, [H, \^ 
or quantum lithography [4l|, [H, [4^ . For advanced ap- 
plications, multimode entangled states will be required. 
However, the technique used for preparation of arbitrary 
two-mode N-photon states [31] cannot be immediately 
extended to multimode case. 

In this paper, we show that we can exploit experimen- 
tally accessible two-photon entangled states generated in 
the process of spontaneous parametric down-conversion 
to extend the group of states that can be prepared us- 
ing only passive linear optics and single-photon detectors. 
We propose a feasible protocol for generation of an ar- 
bitrary quantum state of two photons in three optical 
spatial modes 



iV-target) = a|200) + /3|020) + 7|002) 
-|-(5|110) + e|101) +ry|011) 



(1) 



where the numbers in each bracket denote the number of 
photons in the first, second and third mode, respectively. 
To prepare an arbitrary state ([1]) we use a multiphoton 
interferometer, which is schematically shown in Fig. [TJ A 
two-mode two-photon state can be prepared in this way 
from a separable state |110) as the input. This approach 
however fails if the number of modes becomes higher. We 
will show that using an entangled input state permits 
us to create arbitrary states of the form ([1]). We will 
then generalize our procedure to conditional generation 
of arbitrary multi-mode entangled two-photon state. The 
scheme is cheap in terms of required resources because it 
involves only two photons and vacuum ancillae which is 
a great practical advantage. 

The rest of the present paper is organized as follows. 
In Sec. II wc prove that separable input two-photon state 
is insufficient for generation of arbitrary multimode two- 
photon state via interference and post-selection. In Sec. 
Ill we describe the scheme for generation of arbitrary 
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FIG. 1: General multiphoton interferometer transforming a 
three-mode input state \tpin) into IV'out) using only vacuum 
ancillae. We post-select cases when both photons are present 
in the first three output modes. 



three-mode states ((T]) . Results of numerical simulations 
of generation of various states are presented in Sec. IV. In 
Sec. V we generalize the state preparation procedure to 
arbitrary number of modes. A simple experimentally fea- 
sible scheme tailored for generation of a specific class of 
four-mode two-photon states, the so-called KLM states, 
is discussed in Sec. VI. Finally, Sec. VII contains a brief 
summary of the main results and conclusions. 



II. MOTIVATION 

First of all we prove that a separable input state is not 
sufficient for preparation of an arbitrary state of the form 
(HI) by the interferometric method sketched in Fig.[TJ We 
express the input state |110) as alajlvac), where a| de- 
notes creation operator for the i^^ mode and |vac) stands 
for the vacuum state. General interferometer performs 
linear transformation that can be written as 



(2) 



where Uij € C are elements of a unitary matrix U. For 
future reference, we note that the interference of modes 
a and 6 on a beam splitter BS with amplitude transmit- 
tance t and reflectance r ^ \/l — t"^ is governed by 



= ta. 
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and a phase shifter PS shifts the phase of a single mode 



be transformed as follows, 



according to a^^ 



1110) 



The separable state |110) will 
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FIG. 2: Initialization phase of the protocol. Pair of pho- 
tons obtained by spontaneous parametric down-conversion in 
a nonlinear crystal BBO is sent to polarizing beam splitters 
PBSl and PBS2 which transmit horizontally polarized pho- 
tons and reflect vertically polarized photons. Half wave plates 
HWPl and HWP2 are inserted to modes 2 and 3 to align the 
polarization of all beams to vertical. The photons then inter- 
fere on two ordinary balanced beam splitters BSl and BS2. 



Conditional projection of all output ancillae onto vacuum 
then gives 



1110) 




vac) 



(5) 



k=l 



It is easy to verify that the class of achievable target 
states ([5]) is limited. For example it can be shown that 
the state 



|110) + |101) + |011) = (at^t + fltct + 5tct)| 



vac) 



(6) 



cannot be expressed in the factorized form ([5]) and is 
therefore not obtainable from separable state like |110) 
by multiphoton interference using only vacuum ancillae. 
Addition of single photon ancillae would make the prepa- 
ration possible, but also experimentally more difficult 
0. 

Even with one fixed entangled state as the input, one 
can not prepare arbitrary state ([1]) using only determinis- 
tic unitary transformations. To prove this, let us consider 
two photons in n optical modes. Such quantum state is 
then determined by 



2i?^ - 2 = 72^ 



(7) 



real parameters where i?J! — ("^^ "'^) denotes the num- 
ber of combinations with repetitions of n elements in k 
classes. The term —2 corrects the number of parame- 
ters due to normalization and irrelevant overall phase of 
the state. On the other hand, n-dimensional special uni- 
tary transformation represented hy nx n complex matrix 
contains 



(8) 



real parameters, where = (^) denotes the number 
of combinations without repetitions of n elements in k 
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FIG. 3: Preparation phase. The linear interferometer involves three unitary operations UOll, UlOl and U(BS3) composed 
of beam splitters BSj and phase shifters PSk. The preparation also requires two attenuation filters Fl and F2 which are 
implemented using unbalanced beam splitters and vacuum ancillae. Successful filtration can be verified by monitoring the 
output ancilla modes with single-photon detectors. In the design of the scheme we imagine reverse propagation of the state 
backwards from the output to the input. 



classes. The factor ~n in Eq. ([8]) accounts for the reduc- 
tion of the number of free parameters due to normaliza- 
tion condition of each row of the matrix. The term — 2C^ 
represents correction due to condition of orthogonality of 
rows and the factor —1 occurs since the determinant of a 
matrix belonging to the SU(n) group is fixed and equal 
to 1. It is evident that for any n > 1 the number of pa- 
rameters specifying the quantum state is larger than the 
number of degrees of freedom allowed by unitary trans- 
formation. This results in the need of probabilistic non- 
unitary filtering operations in a fully general preparation 
scheme. 



III. PREPARATION SCHEME 

Here we show that by exploiting the two-photon entan- 
gled states generated by spontaneous parametric down- 
conversion as a resource, we can conditionally generate 
arbitrary state of the form ([T]) with only linear optics. 
The protocol can be divided in two phases: initialization 
phase and preparation phase. We begin by describing 
the first one. 



A. Initialization phase 

We use spontaneous parametric down-conversion to 
prepare an entangled pair of photons (see Fig. [2]). The 
initial state of these photons can be written as 



1 



IV') = ^ (l^^v^) 



\VH))., 



(9) 



where H and V denote state of single photon polarized 
horizontally or vertically, respectively. Using two polariz- 



ing beam splitters that reflect vertical and transmit hor- 
izontal polarization we are able to obtain the four-mode 
state 



1 



1^) = ^ (11001) + |0110», 



(10) 



where the Fock state basis is used. Four numbers in each 
bracket denote the number of photons found in spatial 
modes 1 through 4, respectively. Half-wave plates HWPl 
and HWP2 unify polarization states in all four modes. 
Then the Hong-Ou-Mandel interference [1^ takes place 
on balanced beam splitters BSl and BS2 and the state 
of the photons changes to 

1^) = i (|2000) + 10200) -t- 10020) -I- |0002)) . (11) 

If we want to generate a three-mode state, then we simply 
omit the last mode. Since in our protocol we post-select 
cases when two photons are present at the output of the 
preparation device, we can assume the effective initial 
two-photon state in the form 



1 



|V') = ^(|200) + 1020) + 1002)), 



(12) 



which will be used in the second phase of the protocol. 



B. Preparation phase 

In order to prepare an arbitrary state H]) the state 
is fed into a multiport linear interferometer schematically 
illustrated in Fig. [31 The device is composed of passive 
linear optical elements such as beam splitters and phase 
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shifters and involves filters FX and F2 where the signal 
beams are attenuated by mixing them with auxiliary vac- 
uum beams on beam splitters. The filtering operation 
is successful if the photons do not leak into the output 
auxiliary ports. This can be, in principle, verified by 
means of single-photon detectors placed on the auxiliary 
outputs, as shown in Fig. [31 However, such verification 
would require perfect detectors with unit efficiency. In- 
stead, we can utilize a simpler verification strategy rely- 
ing on postselection. To prove that the scheme is work- 
ing, one can set detectors on every output mode and look 
for double coincidences which herald successful prepara- 
tion of the state. Quantum state tomography [46| could 
be used to completely characterize the generated state. 
In this case the target state is destroyed by the measure- 
ment. If the prepared state serves as an input for another 
scheme, one can post-select successful preparations at the 
output of such scheme. 

The parameters of the interferometer elements can be 
determined by "reverse engineering" of the target state. 
We let the target state propagate through the scheme in 
reverse and we successively get rid of all terms but the 
terms that are present in the state (fT^ . The reversal is 
easily done for beam splitters (BSl, BS2 and BS3) and 
phase shifters (PSl, PS2 and PS3) because they are rep- 
resented by unitary transformations. Filters Fl and F2 
are also (probabilistically) reversible. Once the reverse 
procedure finds correct parameters for all elements, the 
scheme is designed to generate the target state ([T]) from 
the state ([T^ . 

In the Heisenberg picture, passive linear transforma- 
tions can be described as linear unitary transformations 
of the creation operators, c.f. Eq. Every three-mode 
transformation can be decomposed into a sequence of 
two-mode and single-mode transformations correspond- 
ing to beam splitters and phase shifters [13]. A fully 
general interferometric scheme can thus be constructed 
from these basic building blocks. Such a device would 
provide a large number of degrees of freedom that can 
be exploited to optimize the success probability of the 
state preparation but at the same time its experimental 
realization would be quite difficult. We have therefore 
chosen some sort of compromise. Our scheme includes 
all necessary elements needed for the preparation of ar- 
bitrary state H]) and one additional beam splitter BS3 
and phase shifter PS4 as optional optimization elements. 
These components can be removed and our scheme thus 
simplified at the expense of a reduced probability of suc- 
cessful state preparation. 

In what follows we imagine that the state propagates 
backwards through the scheme from the right to the left. 
We successively eliminate all terms where two photons 
are present in two different spatial modes and end up 
with the state as a result. First of all we get rid 
of the cross term |101) with complex amplitude e by de- 
structive interference. To this end, we employ a unitary 
transformation UlOl represented by a beam splitter BS2 
with amplitude transmittance T2 accompanied by phase 



shifter PS3 imposing phase shift cjj^. By a simple calcu- 
lation, one finds that setting c/jt, — — argi5 -I- arge and 
T2 = -^^=^, where u= |||, the state |101) is completely 
eliminated by destructive interference. 

The success rate of the protocol can be increased by 
the optimization beam splitter BS3 mixing modes 1 and 
2. Amplitude transmissivity of BS3 should be chosen 
such that in the reverse propagation after transformation 
UlOl the amplitude of the state |110) is minimized. The 
optimal value of T3 can be found numerically by a simple 
computer algorithm that maximizes the overall probabil- 
ity of success of the protocol. The reverse propagated 
state impinging on the filter F2 from the right can be 
written as 

\^') = a'|200) -I- P'\020) + 7'|002) + +6'\110) + ry'|011). 

(13) 

If the term |110) is not completely eliminated and 5' ^ 
we are forced to employ a filter F2 in modes 1 and 2. 
Let and 6^ denote the creation operators for modes 1 
and 2, respectively. The filtration by F2 can be expressed 
as 



'^/2 



92 



(14) 



where is set so that the matrix can make part of a 
bigger unitary transformation on three modes. Without 
loss of any generality we can assume that b' j a' is real and 
positive because the phase shift between the two ampli- 
tudes is compensated by the phase shifter PS2 if we set 
02 — arga — arge. The amplitude transmittance of the 
beam splitters forming the filter F2 is equal to q^. The 
correct mixing of modes is achieved if the product of the 
amplitude reflectances \ — q^oi the two beam splitters is 
equal to —b'qil\/2oL. This yields a quadratic equation 
for qi whose solution reads 



92 



V2o 



a' 



(15) 



The filtering requires an ancilla vacuum mode and the 
transformation p4)) is effectively implemented if no pho- 
ton leaks into the output ancilla. Note that Eq. (fH)) 
does not preserve canonical commutation relations be- 
cause the operator of ancilla mode is omitted. This sim- 
plified mathematical description of the filter is neverthe- 
less correct for our purposes. 

After filter F2 the state simplifies to 

IV") = a"|200) + /?"|020> + 7"|002) + 77"|011). (16) 

To ensure correct state preparation, one has to compen- 
sate filtering of modes 1 and 2 by additional attenua- 
tion of the third mode by a beam splitter with amplitude 
transmissivity q^. The probability of success Pp2 of the 
filter F2 can be lower bounded as follows. The success- 
ful filtration during the state preparation procedure can 



5 



be described by the transformation = M|V'"): where 
the non- unitary operator M depends on g2- We have 
Pf2 = (V'"|MtM|?/'")/(V'"IV'") from which we obtain the 
inequahty Pf2 ^ Wmin where m„iin is the lowest eigen- 
vahie of the matrix M'I'M. A straightforward calculation 
yields 



(1 



2 
92 



4f] 



4- 



(17) 



It can be shown that m„iin is a monotonically increasing 
function of (72 which suggests that in order to maximize 
the success rate of the protocol we should maximize 52- 
It follows from Eq. that monotonically decreases 
as the ratio 5' /a' increases. It is therefore desirable to 
minimize (5'and maximize a' which is precisely the pur- 
pose of the beam splitter BS3. Let us now derive a lower 
bound on g|. The maximum of absolute values of the 
amphtudes specifying the state H]) is greater or equal to 
If the maximum amplitude corresponds to one of 
the states with two photons in a single mode then by 
simple relabeling of the modes we get a' > l/\/6- If 
this amplitude corresponds to mode with two photons in 
two different modes, say |110) after possible re-labeling, 
then we may use balanced beam splitter BS3 together 
with proper phase shifts to reach amplitude a' > 1/^/12. 
Since the unitary operations U(BS3) and UlOl are deter- 
ministic, the state (1131) is normalized and la 



which implies \S'\ < Taking everything together. 



we find that |(5'/a'| < VH can always be satisfied so that 
g| > (15- /4 w 0.136 holds. If we plug in this lower 
bound into Eq. ^ we finally get Pf2 > 0.52%. 

In the next step of the protocol, the term 77"|011) in 
the state ([T6)) is removed by a destructive interference be- 
tween second and third modes. This is accomplished by 
a unitary transformation UOll which consists of a phase 
shifter PSl imposing phase shift 0i on the second mode 
and a beam splitter BSl with transmittance ti = cos?? 
mixing the second and third mode. The term ry"|011) is 
eliminated provided that 



tanc^ ~ 



tan(2z?) 



Im[(7" - /3")/r?"] 
Rc[(7" + /3")/r]"] 
V2rj" 



After the whole procedure we obtain the state 
\^"') = a"'|200) + /3"'|020) + 7"'|002), 



(18) 



(19) 



which can be transformed by filters Fl into the state 
(fT2)) . We denote by Ppi the overall success probability 
of filters Fl. In order to obtain the state from the 
initial state (|12p we in fact need to apply attenuation 
filters to only two modes. With probability i the two 
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FIG. 4: Success probability Ps of the protocol for states (|2ip 
is plotted as a function of the parameter w € [0; ^] which 
specifies the "difi'erence" of the target state from the initial 
state dnj. 
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FIG. 5; Success probability Ps of the protocol for states (p2l 
is plotted as a function of the phase shift g £ [0, 2tt]. 



photons are present in a mode which is not attenuated 
by Fl and, consequently, it holds that Ppi > i. 

Parameters of all components are now determined. By 
proper inversion of the filters we obtain configuration 
which prepares the target state ([T|) from the state (|12p . 
The total success probability for the whole protocol is 
defined as 



-Ps = (?Aout|V'out), 



(20) 



where IV'out) is un- normalized output state obtained 
from the normalized input state (HI]). It holds that 
Ps = PfiPf2 > 0.17%. 



IV. NUMERICAL SIMULATIONS 

To verify the functionality of our protocol, we have per- 
formed extensive numerical simulations of preparation of 
various states. An example of the results is given in Fig. 2] 
which displays the probability of successful preparation 
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FIG. 6: Success probability histogram. For a given randomly 
generated set of 10^ states, the minimum success probability 
of the scheme reads 0.69%. The average success probability 
is around 26.0%. 

for target states 

cx cosw(|200) + |020) + |002)) 

+ sini«(|110> + |101> + |011>), (21) 

where w S [0,-|]. We can see that as w increases the 
target state is "shifting off" the state (IT^ which is the 
product of the initiahzation phase. Even though the 
probabihty decreases, it is stiU weU above 30%. It can 
be shown that for this particular class of states the best 
performance of the protocol is obtained by using a bal- 
anced optimizing beam splitter BS3. This probability Ps 
is however very sensitive to exact splitting ratio of the 
beam splitter when w ~ O.Stt. It this case it is therefore 
more convenient for the experimental realization to use 
an unbalanced beam splitter with splitting ratio that cor- 
responds to another local maximum of the success prob- 
ability. Even though this choice gives a slightly lower 
success probability, it makes the scheme more robust and 
tolerant to imprecision in the splitting ratio. 

In another numerical simulation we have studied 
preparation probability Ps for a more asymmetric single- 
parametric class of states 

IV2) oc |200) + |101> + |002) 

-f-e^'dllO) + |020) + |011)), (22) 

where g E [0, 2tt]. The result (see Fig.[5]) shows a symetric 
Gauss-like function realising maximum of 48% at g = tt 
and asymptotically aproaching minimum of 34% for g = 
0,27r. 

Finally, we present in Fig. [5] a histogram of the suc- 
cess probability Ps of the preparation protocol. This his- 
togram was obtained by evaluating Ps for 10^ randomly 
generated states. We can see that most of the states can 
be prepared with probability larger than 15%. However, 



there is a small class of states that are hard to generate. 
The lowest success probability obtained by the numerical 
calculations reads 0.69% which qualitatively agrees with 
the lower bound 0.17% derived in the preceding section. 

V. GENERALIZATION TO AN ARBITRARY 
NUMBER OF MODES 

So far we have explained how linear optics can be used 
to prepare an arbitrary state of two photons in three spa- 
tial modes. We now generalize our procedure to an arbi- 
trary number of modes. Figure [7] schematically depicts 
the scheme for preparation of an arbitrary two-photon 
state of the form 

N N 
l^target) 12:^)+ E M'^M- (23) 

Here N denotes the number of modes, \2j) represents 
state with two photons in jth mode and vacuum in all 
other modes and |ljlfc) indicates state with one photon 
in modes j and k an no photons in the other modes. 

We begin by employing non-degenerate spontaneous 
parametric down-conversion to produce correlated state 
of two photons in fi = N/2 pairs of modes ja and jb, 

1 

IV-SPDC) = -/=$]|l,a)|l,fc). (24) 

i=i 

The modes ja and jh can be for instance different direc- 
tions on the cone of the emission from the crystal [isj . 
Alternatively, one can pump the non-linear crystal with 
a sequence of d ultra-short pulses and then the modes 
ja and jh correspond to jth time-bin [49[. Yet another 
option is to simultaneously pump d nonlinear crystals. 
After down-conversion, the pairs of modes ja and jh are 
combined on d balanced beam splitters and the Hong- 
Ou-Mandel interference produces the state 

1 ^ 

In the second part of the protocol, the A^-mode two- 
photon state (|25p is injected into the interferometer 
shown in Fig. [T] The parameters of the interferometer 
can be again determined by reverse engineering where 
we let the target state propagate backwards through the 
scheme and at each step we eliminate certain component 
of the state by destructive interference. In this way we 
finally obtain the state ((25|l . We now describe the cy- 
cles of the protocol in more detail. In the first cycle, we 
get rid of all states with exactly one photon in the last 
mode. First we apply unitary transformation to the first 
and second mode. This transformation is equivalent to 
the transformation UlOl from Fig. [3] and eliminates the 
term |100...001). Next the second unitary transformation 
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FIG. 7: Scheme for preparation of generic A^-mode two-photon state. As explained in the text, the scheme works in cycles. 
The scheme is composed of unitary transformations U that consist of a beam splitter and two phase shifters before and after 
the beam splitter, filters F equivalent to filter F2 from the scheme in Fig. [S] and the optimizing beam splitters O. 



acting on the second and third mode eliminates the term 
|010...001). We continue in this fashion with only one 
exception. We do not use unitary transformation to mix 
the (A^-l)th and iVth mode. Instead, we need to use a 
filter F similar to the filter F2 shown in Fig. [3l accompa- 
nied by attenuation of the first N ~ 2 modes (not shown 
in the figure). This filter removes the term |00...11). At 
this stage, all terms with exactly one photon in the last 
(iVth) mode are eliminated and there is either zero or 
two photons in the last mode. For subsequent considera- 
tions it can be therefore neglected and we have effectively 
reduced the A^-mode problem to (A^-I)-mode problem. 

In the second cycle we apply the same strategy to 
eliminate all cross terms with one photon in the mode 
A^— 1. Repeating this procedure N—3 times we finally get 
back to the three-mode situation which we have already 
solved. The success probability of the protocol can be 
improved by employing optimizing beam splitters before 
each cycle. These beam splitters play the same role as 
BS3 in Fig. [3] and they should combine the modes where 
the filter will be applied in a given cycle. The number 
of required two-mode operations that have to be imple- 
mented scales as ^N{N — 1) with N being the number 
of modes. 



VI. PREPARATION OF TWO-PHOTON KLM 
STATES 

The general scheme presented in the previous section 
can prepare any two-photon quantum state in an ar- 
bitrary number of modes. The experimental difficulty 
however increases considerably as the number of modes 
grows. If one is interested in preparation of specific class 
of quantum states then one could attempt to design a 
less general but simpler and experimentally more feasi- 
ble scheme for this purpose. As an example of an impor- 
tant specific class of states we consider here the so-called 
two-photon four-mode KLM states, 

IV'KLm) = a|f fOO) -I- /?|0f fO) -I- a|00f f>. (26) 




PBS 



FIG. 8: Scheme for preparation of the two-photon four-mode 
KLM states. 



These states were introduced by Knill, Laflamme, and 
Milburn in the context of quantum computing with lin- 
ear optics [so*! . They can be used for teleportation-based 
implementation of quantum gates and the teleportation 
fidelity can be optimized by appropriate tuning of a/P 
ratio [5l| . In view of these applications it is clearly im- 
portant to be able to prepare such states in the simplest 
possible way. 

The scheme for preparation of the two-photon KLM 
states is shown in Fig. [8l A strong coherent laser 
beam pumps nonlinear crystal BBO cut for Typc-II 
phase matching where a maximally entangled two-photon 
triplet state = -^{\HiV2) + \V1H2) is generated in 
the process of spontaneous parametric down-conversion 
(indexes 1 and 2 denote the spatial modes). The first 
spatial mode is then subjected to polarization sensitive 
filtering using a beam splitter F which entirely transmits 
vertical polarization and has amplitude transmissivity r 
for horizontal polarization. In practice, such filtering can 
be accomplished, e.g., by means of tilted glass plates [s^ ]. 
The two spatial modes are then combined on an ordinary 
(polarization in-sensitive) unbalanced beam splitter BS. 
We parametrize the transmittance t and reflectance r of 
BS by an angle and we have t = cosi? and r = sini?. 
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FIG. 9: The success probability of the protocol (full line), the 
probability of finding the state |0110) in the KLM state (dot- 
ted line) and the amplitude transmissivity of the beam splitter 
BS (dashed line) are plotted as functions of the parameter r. 

In order to generate the desired state ((26)) we have to set 

iJ^arctanVr. (27) 

The polarization sensitive phase shifter PS adds a phase 
shift of TT to the vertically polarized photon in the first 
output spatial mode. As a result of this entire procedure 
we obtain the state 

1^) - a\HiVi) + (3\Hi^2) + a\H2V2), (28) 

which can easily be split by two polarizing beam splitters 
to obtain the state ([26| with parameters 

Since the scheme contains a filter F the state is prepared 
only probabilistically and the success probability is equal 



to Ps = ^(1 + T^)- Figure [9] illustrates the trade-off 
between success probability of the whole protocol and the 
probability of finding |0110) in the generated state. 
The teleportation protocol proposed by KLM requires 
the state (pB)) with a = (3 and for this case we find Ps = 
0.57. 



VII. CONCLUSIONS 

In this paper we have demonstrated how entangled 
pairs of photons generated in the process of spontaneous 
parametric down-conversion can be used for preparation 
of an arbitrary multimode two-photon state. This prepa- 
ration requires only linear optical elements (beam split- 
ters, wave plates, phase shifters) whose parameters can 
be determined analytically by reverse engineering proce- 
dure whereby the target state is propagated backwards 
through the scheme. We have also suggested a simple 
and experimentally feasible setup specifically tailored for 
preparation of two-photon KLM states which are cru- 
cial for linear optics quantum computing. The proposed 
schemes rely only on existing and well mastered exper- 
imental techniques and can be therefore implemented 
with current technology. Our findings thus pave the 
way towards generation of complex multimode entangled 
states of light required for advanced quantum informa- 
tion processing, ultraprecise measurements or quantum 
lithography. 
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